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find manifestations of non-linear screening and charging effects around integer filling 
factors, consistent with recent imaging experiments. Our calculations exhibit g-factor 
enhancement as well as strong overscreening in the centre of the Landau bands. Even 
though the critical behaviour appears mostly unaffected by interactions, important 
implications for the phase diagram arise. Our results are in very good agreement 
with the experimental findings and strongly support the relevance of electron-electron 
interactions for understanding integer quantum Hall physics. 
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1. Introduction 

The integer quantum Hall effect (IQHE) — observed in two-dimensional electron 
systems (2DES) subject to a strong perpendicular magnetic field [1] — has been 
well explained based on single-particle arguments [2-11]. Crossing the centre of each 
disorder- broadened Landau level, the Hall conductivity axy jumps by e^/h and the 
longitudinal conductivity a^^ is finite. The accompanying localization-delocalization 
transition is governed by a critical exponent i) = 2.34 ±0.04 [12,13], a universal quantity 
independent of microscopic details of electron motion, the disorder realization, and thus 
the chosen material. In fact, it is precisely this astonishing resilience of the quantum 
Hall (QH) effect that makes it an ideal metrological standard [14]. 

However, recent experiments on mesoscopic MOSFET devices have questioned the 
validity of such a simple single-particle picture. Measurements of the Hall conductance 
as a function of magnetic field B and gate voltage [15] exhibited regular patterns along 
integer fiUing factors. It was argued that these patterns should be attributed to Coulomb 
blockade effects. Similar patterns have been found recently also in measurements of the 
electronic compressibility k as a function of B and electron density Ue in the IQHE [16] 
as well as the fractional quantum Hall effect (FQHE) [17]. From these measurements it 
turns out that deep in the localized regime between two Landau levels, stripes of constant 
width with particularly small k can be identified. These stripes consist of a collection 
of small-K lines, identifiable with localized states, and their number is independent 
of B. This is inconsistent with a single-particle picture where one expects a fan- 
diagram of lines emanating from (0, 0) in the {B, ne)-plane. The authors argue that their 
results may be explained qualitatively by non-linear screening of the impurity charge 
density at the Landau level band edges. Clearly such screening effects — explained 
within a Thomas-Fermi approach [16] — are beyond a simple non-interacting theory. 
This immediately raises a question on the status of the aforementioned universality 
of the QH transition [15, 18, 19], which was obtained largely within a single-particle 
approach [9,12,20-25]. Further qualitative support for the screening arguments where 
recently given in [26,27], where a Hartree-Fock (HE) approach lead to a non-fan structure 
for K in the {B,ne) diagram; a further Thomas- Fermi- type investigation showed that 
charging lines with the desired slopes can indeed be found within a simple model for 
the impurity potentials [28]. 

In the present paper, we quantitatively investigate the effects of Coulomb 
interactions on the compressibility in the (i?,ne)-plane within a HE approach. HF 
accounts for Thomas-Fermi screening effects while at the same time leading to a 
critical exponent u whose value is consistent with the results of the non-interacting 
approaches [29,30]. We find that the observed charging lines in the compressibility can 
be well reproduced and that the width of each group of lines is well estimated by a 
force balance argument. Thus we can quantitatively explain the lines and stripes in the 
compressibility as a function of (B^rie). We note that these results fully support the 
qualitative picture proposed in [16,26-28]. 
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2. Formulation of the QH model in Landau basis 

In order to model a high-mobility hctcrostructure in the QH regime, we consider a 2DES 
in the {x,y)-p\ane subject to a perpendicular magnetic field B — Be^. The system can 
be described by a Hamiltonian of the form 

^2DES = /i^ + T^c = + 2 ^ 

where cr = ±1 is a spin degree of freedom, is a smooth random potential modeling the 
effect of the electron-impurity interaction, Vq represents the electron-electron interaction 
term and m*, g*, and //b are the effective electron mass, g^-factor, and Bohr magneton, 
respectively. In order to avoid edge effects wc impose a torus geometry of size L x L 
onto the system [31]. The electron-impurity interaction is modeled by an electrostatic 
potential due to a remote impurity density separated from the plane of the 2DES 
by a spacer-layer of thickness d, as found for instance in modulation-doped GaAs- 
GaAlAs hetero junctions. Within the plane of the 2DES, this creates a random, spatially 
correlated potential with a typical length scale d. We use A^i Gaussian-type " impurities" , 
randomly distributed at r^, with random strengths Ws G [—W, W], and a fixed width d, 
such that Vi{f^ = Efii (ws/^d^) exp[-(f - fsf/d^] = J2^VM exp(i? • ^ with 

VM-f:^.e.Jy^-^,^.r^, (2) 

s=l ^ ^ 

where q^^y = 2TTj/L and j = —N^, —N^ — 1, ... , A*",^. The areal density of impurities 
therefore is given by ni = Ni/L"^. The limit (i ^ yields a potential of 5- 
type that would be more adequate for modeling low-mobility structures. We will 
highlight a few differences between the two cases in the following sections. The 
electron-electron interaction potential has the form Vc{r,r') — 7e^/47reeo|r — r*| = 
Hq^ciq) exp [iq- (f - r')] , with 

The parameter 7 will allow us to continually adjust the interaction strength; 7=1 
corresponds to the bare Coulomb interaction. Choosing the vector potential in Landau 

— * 

gauge, A = BxCy, the kinetic part of the Hamiltonian is diagonal in the Landau 
functions [32] 

1 r, ix-kiiY' 

iky — 



^n,k[r = , =exp 



C I 



Hn [ ) , (4) 



211 

where n labels the Landau level index, k — 27rj/ L with j — 0, . . . , A^ — 1 labels the 

momentum, Hn{x) is the nth Hcrmite polynomial, and 1^ = \Jh/eB the magnetic length. 
These functions arc extended and L-periodic in y-direction and localized in x direction. 
For the system's many-body state, |$), wc use the usual ansatz [33,34] of an anti- 
symmetrized product of single particle wave-functions 4>a{f) (Slater determinant), which 
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we choose as a linear combination of Landau states 

= E E CZkXn,k{rl, (5) 

n=0 fc=0 

with A^LL being the number of Landau levels and the periodic Landau functions 

oo 

Xn,fe(r) = (f^nk) = E fn,k+jL/pS'^, (6) 

j = -OD 

in order to meet the boundary conditions. The number of flux quanta piercing the 
2DES of size L x L is given by = L^/27rZ^, yielding a total number of M = A^ll^^ 
states per spin direction. The filling of the system is characterized by the filling factor 
u — Ne/N^, with A^e being the number of electrons in the system and areal density 
Tie — Ne/L"^. The total Landau level density is given by no — eB/h. 



3. Hartree-Fock equation in the Landau basis and its numerical solution 

We arc left with finding the correct expansion coefficients C"'^ [33-36]. The Hamiltonian 
is represented in matrix form using the periodic Landau states \nk) and we have 

ii:,k;n',k' = {nk\H^^^,\n'k') 

with the cyclotron energy hwc = hcB/m*. The disorder matrix elements are given by 
^n,k;n',k' = J^q^ii.'l) ^n,k;n' ,k' {(f) , whcrc mixing of Laudau levels is included. The exphcit 
form of the plane wave matrix elements Sn,k;n',k'{Q) = {nk] exp(ig • r)\n'k') is computed 
in the Appendix. The elements of the Fock matrix F are 

xpo- \ ^ \ / ^m,l;m',l' _ r ^m,l;n' ,k'\ j^a' (q\ 

^ n,k;n',k' — ^ \^ n,k;n' ,k' "<^,<^'^n,k;m',l' J ^m,l;m',l' ' \°) 

a' l,Tn,l',Tn' 

where the first term is the Hartree and the second the Fock contribution. The 
bielectronic integrals G^^p^fc, = J2^^oVc{q)Sn,k;n',k'{q)Sm,i;m',i'{-q) can be further 
simplified as given in the Appendix. A homogeneous, positive background is assumed 
that neutralizes the charge of the electrons and thereby prevents the Coulomb term 
from diverging as |^ — > 0. In fact, this interaction with the background can be shown 
to cancel with the term |g| = in F up to a contribution of the order of due to the 
finite system size [37] . The density matrix is given by 

M 
a=l 

with Tr(D) = A^e and D'^D^ = D'". Here /(e) denotes the Fermi function. The total 
energy E'tot in terms of the above matrices is given as 

Stot = Tr(hD + ^FD) = iE E i^K,k;n',k'+K,k;n',k')^:,k;n',k'-m 

(J n,k;n' ,k' 
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A variational minimization of (^|i?2DEs|^) with respect to the coefficients yields the 
Hartree-Fock-Roothaan equation [38], a self-consistent eigenvalue problem which in 
compact form can be written as 

= CE", (11) 

with C"' — (Cf , . . . , C^) being the matrix of eigenvectors and E^^ = diag(ej, e^) the 
diagonal matrix of the eigenvalues < ef < . . . < [39] . 

Following the aufbau principle [40], the density matrix is constructed starting from 
the energetically lowest lying state up to the Fermi level ep- In our calculations, we 
keep Nc fixed and compute e-p as the energy of the highest occupied state afterwards. 
Since the Fock matrix depends on the density matrix, which in turn depends on 
the full solution of problem, (11) has to be calculated self-consistently which is 
numerically quite challenging. In the first step we use the solution of the non-interacting 
Hamiltonian h'^ = {nk\h"\n'k') as a starting guess for the coefficients C^. From 
this solution, C^^^^ we construct the density and Fock matrices and finally the full 
Hamiltonian. Diagonalization yields an improved solution, C^^^. The process continues 
until convergence of the density matrix has been achieved. In practice, we compute the 
norm of the difference between successive density matrices ||D'^"'+-'^) — D*^")!! < e. Here 
II • II denotes the Hilbcrt-Schmidt norm defined as ||A|| = Tr(AA*)-'^/^. 

However, convergence of this Roothaan algorithm [38] is not always assured. 
Especially for fillings close to zero or an integer value it may run into an oscillating limit 
cycle. Subtracting a small multiple of the density matrix, feD'^, from the Hamiltonian 
H*^ = h*^ -I- F*^ favours already occupied states, thereby suppressing such oscillations. 
The eigenvalues shift by ef — > ef — 6 V ef < ep, but eigenstates are not affected. 
This algorithm is known as the level shifting algorithm [41] and h is the level-shift 
parameter. We choose h ^ 1.76e^/L which is exactly the contribution of the neglected 
Vc{q = 0) term that does not cancel with the background charge due to the finite system 
size. This prevents the occupied and unoccupied states from mixing and stabilizes the 
convergence of the algorithm, although at the cost of slower convergence. A further 
improvement of convergence is achieved by constructing the new density matrix D^""^-^^ 
as a mixture of the updated and the previous one, i.e. D^^^^^ = AD^^^^^ + (1 ~ A)D("), 
with A e [0,1]. We adjust A in each step automatically such that the fastest global 
convergence is guaranteed [40]. 

In each HF step, assembhng the dense Fock matrix F*^ scales as 0{Nlj^N^) and is 
clearly very time-consuming. An improved scheme, even though generally possible, is 
of little advantage since the diagonalization is of similar complexity. For the calculation 
of a particular disorder configuration and magnetic field, a self-consistent run has to be 
made for each of the M possible filling factors. Hence, the complexity of a complete 
HF calculation is of the order 0{2KNl^N^) with K the number of iterations until 
convergence. The dependence on the system size is therefore 0{L^). For system sizes 
of L ~ 300nm, we find K ~ 100 — 1000. In all results present here, convergence of the 
HF scheme is assumed for e < 10~^. 
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4. Chemical potential and electronic compressibility 

The electronic compressibility k = {drie/ dfi) /n1 reflects the ability of the 2DES to 
absorb electrons when changing the chemical potential. With ji — dE^ot/dNe, we find 
dji/drie — L'^{d'^Etot/dN^). Hence, for finite sample calculations, we can obtain k, from 
Etot(Ne) using 

§^^L^ [EtotiN, + 1) - 2Etot(iVe) + ^tot(A^e - 1)] • (12) 

Alternatively, at T = 0, we can compute the change in the chemical potential for Ne 
electrons by noting that the Fermi energy eF(A^e) = A*(A^e)- Thus we immediately have 

^ = L' [e^iN + 1) - e^N)] . (13) 

This turns out to be numerically more stable than (10) and shall be used in the following. 

When the Fermi energy lies in a region of highly localized states, it takes a 
considerable energy to accommodate another electron, and thus the compressibility is 
low. On the other hand, in a region of delocalized states a newly added electron is much 
more easily absorbed and k is high. For a non-interacting system, n is proportional to 
the tunneling density of states (DOS) [42] and expected to exhibit a fan-hke structure 
in the {B, rig) diagram. In particular, the resonances in k need not align with slopes 
equal to integer filling factors [43]. In the interacting case, k, is proportional to the 
thermodynamic density of states (TDOS) and the inverse screening length. We note 
that experimentally the change of the chemical potential is detected when changing the 
back gate voltage and hence the electron density at constant B [16, 17]. 

5. Sample mobility and DOS 

In order to tune our parameters for the electron- impurity interaction to the experiment, 
we estimate the zero field mobility, defined as /Iq — er/m* , with r being the transport 
scattering time [42,44]. For a short-range (5-impurity potential, r is identical to the 
single-particle momentum relaxation time, Tg, which determines the level broadening, 
(r/2)^ = huJc.h/2'iTTs [45]. For long-range potential, however, these two times can be 
very different [46] and knowledge about the level broadening does not necessarily imply 
knowledge about the mobility and vice versa. In fact, for a smooth potential with 
d » Zc we have (r/2)^ = ([Vi(f) — (Vi(r))^]^)-. , which does not depend on B. In that 
case, we can determine the mobility from the transport cross-section calculated in Born 
approximation [32,45]. Since the transport scattering time is momentum dependent, 
we take the low temperature limit, where the relevant scattering time is the one for 
electrons having Fermi momentum kp = {27ineY^'^- With a radially symmetric electron- 
impurity interaction potential (for a single scatterer), u{r), we obtain for the transport 
scattering time 

= "^J^eil-cosiO)] \u{2kFsm{9/2))\\ (14) 
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with the Fourier transform u{k) = f (Pfu{f)exp{—if- k) and u{k) = u{\k\). In case 
of 5-interaction and uniformly distributed strengths, we simply have u{k) — VF/-\/3 
and the cos(^^) term in (14) vanishes. Without it, (14) becomes the expression for Tg, 
which proves the equivalence of r and Tg for short-range potentials. For long-range 
potentials, however, forward (small 9) scattering receives little weight since it hardly 
impairs the electron movement and t/ts S> 1. In order to model a situation comparable 
to the experiments of [16], we use material parameters for GaAs {g* — —0.44 for the 
effective g^-factor, m* — 0.067me for the effective mass) and impurity parameters of 
W/nm^ ~ 4eV with a concentration of rij = 3.2 • lO^^cm-^ (e.g. A^i = 288 for L = 300nm 
or A^i = 392 for L = 350nm). Assuming ^ ni, for the (5-potential this yields a 
mobility of ~ lO^cm^/Vs, whereas for d — 40nm (fti spacer layer thickness) we get 
Ho ~ lO^cm^/Vs — a value which is reasonable for a high mobility sample such as a 
GaAs-GaAlAs heteroj unction. 

We have calculated the DOS in the lowest Landau level by averaging over at 
least 1000 samples in the non- interacting and the interacting system. In Figure 1 
we show results for Gaussian and 5-impurities for three system sizes L — 400, 500, 
and 600nm for the non-interacting case as well as the interacting case. Within the 
accuracy of the calculation, we find the DOS to be independent of the system size, 
irrespective of interactions, as of course expected for our static disorder model [47]. 
Furthermore, in the interacting case, we find a strong suppression of the DOS at the 
Fermi level. The formation of this Coulomb gap and its non-criticality has been studied 
previously [29,30,37,48]. In case of Gaussian-type impurities we also observe a strong 
reduction of the band broadening due to screening of the impurity potential, while the 
bandwidth in systems with (5-type impurities is hardly affected by interactions. 

6. Scaling of the participation ratio 

The participation ratio Pa is defined as the inverse of the variance of the charge density 
in the state a, 

Pa-^L'I c^VlVa(r-)r) ' . (15) 

Large values of correspond to spatially extended states, while low values indicate a 
confined state [49,50]. This is intuitively understood by the fact that the density of an 
extended state varies much less over space than a highly localized one. Thus Pa is a 
measure of the degree of localization [51] and may be computed in our model as 

P<^^ Y2 X] ^(C^n,J*C^^fc'(C^J)*C^^^''S'n,fc;n^fcKg)'5'm,^;m^^'(~g)■(l6) 

n,n' ,m,m' q 
k,k',l,l' 

It has been shown that unscreened HF-interactions do not alter the critical exponent 
jy while rcnormalizing the dynamical scaling exponent to ^ = 1 [29,30,53]. As a check 
to our HF results, we calculate Pa of spinless electrons in the lowest Landau level for 
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Figure 1. DOS at B = 3T for the lowest Landau level in a non-interacting (top 
row, 7 = 0) and a HF- interacting (bottom row, 7 = 0.3 at = 1/2) QH system 
for 3 system sizes. The left column shows results for S-type impurities {d = 0) 
with W/nm^ = 2eV, the right column corresponds to Gaussian- type impurities with 
d = 40nm and W/nm^ = 4eV {W/d'^ = 2.5meV). The results in all cases are averaged 
over at least 1000 samples. Error bars arc less than the symbol sizes. Note the strong 
Coulomb reduction of the DOS at the critical energy {ea = Ccrit) in the interacting 
systems. 

the same samples as in Section 5. The participation ratio is expected to obey the single 
parameter scaling form [54] 

P„ = L^(^)-^n(Li/ne.-6e„t|), (17) 

with the anomalous diffusion coefficient D{2) ?a 1.6 — related to the multifractal 
character of the critical states [55-57] — and the critical exponent u « 2.3 [12]. Figure 
2 shows the scaling function for the non-interacting and an interacting system at filling 
factor u — 1/2. The scaling function collapses reasonably well onto a single curve 
for both, non-interacting and HF-interacting systems. We find -0(2) = 1.62 ± 0.10 as 
typical average over both, HF- and non- interacting systems as shown in Figure 3. This 
demonstrates the irrelevance of interactions and the type of disorder for the multifractal 
dimension of the critical state, in very good agreement with previous results [27,30,57]. 
A similar fit in the tails of Pa is numerically less accurate but still yields estimates for 
i> between 2 and 2.4, compatible with the expected value 2.34 ± 0.04 [12,23]. 
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Figure 2. Scaling functions of the participation ratio Pa &t B = 3T for the non- 
interacting (top row, 7 = 0) and the HF-interacting (bottom row, 7 = 0.3 at v = 1/2) 
systems averaged over at least 1000 samples and using D{2) = 1.62, i> = 2.34. The 
left column shows results for (5-typo impurities (d = 0), the right column corresponds 
to Gaussian-type impurities with d — 40nm. Values for L have been scaled by the 
magnetic length. Fluctuations in the tails are due to a smaller number of data points. 

7. Charging lines and stripes in the electronic compressibility 

Let us first briefly recall the experimental results of [16,17] most relevant to the present 
investigation. The compressibilities in the {B, 'ric)-plane (i) exhibit only little variation 
in regions close to the QH transitions at half-integer fllling factors, but (ii) show a 
strong variation between Landau bands at integer fillings which by virtue of the relation 
Tie — veB/h correspond to lines of constant slope. Furthermore, (iii) these regions of 
strong variation seem to have a width which is B and Landau level index independent 
and (iv) within these stripes, thin lines of equal slope jeB/h, j — 0,1,... can be 
identified. In what follows, we have calculated the electronic compressibility as outhned 
in Section 4 in the lowest two, spin-split Landau levels for a sample of size L — 300nm 
at magnetic fields between B — 0.2T and B — 6T. 
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Figure 3. Power-law fit of the system size dependence of Pq, at B = 3T according to 
(17) for non- interacting and HF- interacting systems around = Ccrit- ° and o denote 
S-type impurities {d = 0) whereas * and + show results for Gaussian-type impurities 
with d = 40nm. The error bars for the data points are smaller than the symbol sizes. 
The D{2) values are 1.57 ±0.2 (7 = 0, d = 0), 1.6 ±0.1 (7 = 0, d = 40nm), 1.7 ±0.1 
(7 = 0.3, d = 0), 1.62 ± 0.05 (7 = 0.3, d = 40nm). Values for L are scaled by the 
magnetic length. 



7.1. Compressibilities for the non-interacting system 

Figure 4 shows our results for the inverse compressibihty k,~^ of non-interacting (7 = 0) 
electrons in the two lowest orbital Landau levels, including spin. Darker areas in the plot 
reflect states of higher compressibility, hence a more delocalized regime. Lighter areas 
are more strongly locahzed states. Due to the weak Zeemann splitting, we do not observe 
the two spin bands separately. Rather, both bands remain nearly degenerate and lie 
almost on top of each other. Hence, we only find a single, very strongly incompressible 
region between the first and the second orbital Landau level eik u — 2. This broad hne 
is due to the band gap and the highly localized states at the band edges. Other less 
pronounced lines are visible along different filling factors, seemingly mostly emanating 
from (0,0). Some lines even appear to have a varying slope as shown in the inset of 
Figure 4. We interpret these features as the aforementioned fingerprints of scattering 
resonances in the disorder potential [43] which do not necessarily need to align with 
constant filling factors. Moreover, we clearly observe an increasing number of those 
lines with increasing magnetic field. At u — and 4, the compressibility is again low. 
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Figure 4. Inverse electronic compressibility for a non-interacting system of size 
L — 300nm with disorder strength W/(P = 2.5meV in the (i?, ne)-plane. The color 
scale spans two standard deviations around the average of k~^. The inset shows more 
detailed results for the region marked by a black rectangle. 

7.2. Effects of the Hartree-Fock interaction 

We next include interaction with 7 = 0.3. This is not yet the full 7 = 1 Coulomb 
term, but the results are numerically more stable while at the same time not being 
dramatically different from 7 = 1. Furthermore, 7 < 1 is essentially equivalent to 
increased disorder with the full Coulomb interaction present. 

Figures 5, 6, and 7 show results in the (-B, ne)-plane for an interacting system of 
size L = 300nm with disorder strengths W/d"^ = 1.25, 2.5, and 3.75meV, respectively, at 
fixed impurity range d = 40nm. We observe that the exchange interaction results in an 
effective ^^-factor substantially enhanced from its bare value [58-60], leading to a clear 
separation of the two spin bands. This yields two additional strongly incompressible 
stripes aX u = 1 and z/ = 3, indicated by particularly high values. Quite different 
from the non-interacting case, we find that most of the incompressible lines form groups 
which align parallel in the (5, ?T,e)-plane along integer filling factors. Above a certain 
minimal magnetic field, the width of these groups — the number of the lines — is 
independent of the magnetic field and Landau level, forming incompressible stripes of 
constant width around integer filling factors. Overall, this behaviour is strikingly similar 
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Figure 5. Inverse electronic compressibility for a HF-interacting system of size 
L — 300nm with disorder strength W/cP = 1.25meV in the (i?, ne)-plane. The dotted 
lines show estimates based on a perfect screening condition (see text for details). The 
color scale spans two standard deviations around the average of k^^. The inset shows 
more detailed results for the region marked by a black rectangle. 

to the effects observed in the experiments of [16, 17]. Outside the stripes, there is 
hardly any feature in the compressibihties except directly at the QH plateau-to-plateau 
transitions at half-integer fillings where a small increase in compressibility is discernible. 
In these areas between incompressible stripes, the inverse compressibility tends to have 
a very low or even negative value, which relates to a very high or negative TDOS, 
respectively. This effect has been observed experimentally [61,62] and is a signature 
of the exchange interaction. From the proportionality between compressibility and the 
screening length, we can conclude to observe strong overscreening in the areas of negative 
compressibility. We attribute this to the tendency of the HF-interacting 2DES to form 
a charge density wave [35,63]. Furthermore, when comparing Figures 5, 6, and 7 we find 
that the width of the incompressible stripes increases with increasing disorder strength 
W/d\ 
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Figure 6. Inverse electronic compressibility as in Figure 5 but with W/(P ~ 2.5nieV. 

8. Spatial distribution of electronic density and screening 

The spatial distribution of the total electronic density 

M 

= ^~^Y1 Yl Y'^n,k-y,k'Sn,k;n',kiq)exp{-iqr) (19) 

a n,k,n',k' q 

is readily calculated in our model. It details the screening mechanism by providing 
direct insight into the interplay of disorder and interaction. Let us start at the QH 
transition. Figure 8 depicts the critical charge density at = 1/2 for a non- interacting 
system in units of hq. The contour lines show the impurity potential Vi(r) where the 
critical energy Vi{r) = ep is highlighted by a thick line. The charge density evidently 
behaves according to the semiclassical approximation [23] and follows the equipotential 
lines of Vi(r). For the interacting case, however, we expect Thomas- Fermi screening 
theory to apply [64-67]. The electrostatic potential of the charge density 

m = -^[ dh-^P^ (20) 
^ 47reeo J Ir" - r\ ^ ^ 

and the impurity potential Vi{r) form a screened potential V^cr(r) — -|-e0(r). Here, 
n accounts for the positive background. Since a fiat screened potential is energetically 
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Figure 7. Inverse electronic compressibility as in Figure 5 but with W/(P = 3.75nieV. 

most favourable, one expects to find K,cr(^) = cf for the case of perfect screening. 
However, since fiuctuations of the density, 6n{f) = n{r) — n, are restricted between 
an empty and a full Landau level, i.e. < 6n{r) < no, the screening is not always 
perfect but depends on the fiuctuations in the impurity potential as well as on the 
filling factor [64-66]. The plane can be divided into fully electron or hole depleted, 
insulating regions — where n(r) =0 or n{r) = no, respectively — and metallic regions 
— where n{r) lies in between. Depending on the filling factor, the extent of those regions 
varies. Close to the band edge, insulating regions dominate. Screening is highly non- 
linear and transport virtually impossible. On the other hand, if disorder is weak enough, 
there exists a finite range of filling factors in the centres of each band where metallic 
regions cover most of the sample, percolate and render the whole system metallic. The 
disorder is effectively screened and transport greatly enhanced. In that case, the charge 
density nscr(^ can be obtained by Fourier transforming the screened potential. In 3D, 
this simply leads to the Laplace equation. For 2D, however, one obtains [68] 

nscriq) = -'^^\q\Vi{q) + lynoSg^o, (21) 

where the |g| = term is "perfectly screened" by the positive background and thus does 
not contribute to screening of the impurity potential. In other words, in our model only 
the fiuctuations 6n{r) are essential for screening. Hence, in 2D, a perfectly screening 
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Figure 8. Spatial distribution of non- interacting electron density n{r)/nQ sA B — 4T, 
7 = and v = 1/2 as indicated by the color scale. Solid contour lines show the 
equipotential lines of the Vi{r). The thick solid lines corresponds to ep- 



charge density would obey 

47reeo 



A2DVi(r^ 



(22) 



Clearly, the actual charge density is expected to deviate from riser (^) for several reasons. 
Firstly, the fluctuations of n{r) are restricted as discussed above. Secondly, (22) is valid 
for the Hartree case only. Taking the Fock contribution into account will introduce short 
wavelength fluctuations due to the tendency for crystalization. However, we still expect 
the charge density to follow (22) in the limit of |g| ^0 [69]. Figure 9 shows results for 
the charge density of interacting electrons at z/ = 1/2. Broken lines indicate the regions 
where nscr(r) exceeds the range for 5n{f) either below or above, i.e. areas that cannot 
be screened at all and thus exhibit insulating behaviour. Otherwise, we find the charge 
density to follow nscr(^) very closely. In this regime, the density is well described by (22) 
and the screening is very effective. Metallic regions dominate over insulating ones and 
transport is expected to be good. In Figure 10 we depict cross-sectional plots of n{r) 
and riser (0 for the sample of Figure 9 at a; = lOOnm and three different filling factors, 
demonstrating the discussed effects again very clearly. 
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Figure 9. Spatial distribution of HF-interacting electron density n{r)/no at B — 4T, 
7 = 0.3 and i/ — 1/2 as indicated by the color scale. Contour lines show (21). The 
broken lines indicate unscreenable (insulating) regions. The thick solid line shows 

"scr(^) = ^scr = 0. 

9. Breakdown of linear screening 

Thus far we have shown that our results can qualitatively reproduce the structures 
observed in the {B,ne) plots of the compressibility. We find stripes of constant width 
with very similar characteristics as in the experiments. Furthermore, we show that 
within HF, the impurity potential in the band centre can be quite effectively screened 
by the charge density. Let us now turn our attention to the stripes. The screening of 
the impurity potential is non-linear near the edges of the Landau bands. Most of the 
sample is thus covered by insulating regions where the Landau band is either completely 
depleted or filled. Metallic behaviour is confined to small regions around potential 
extrema, where electron or hole islands are formed. If additional charge is introduced 
into the system, charging effects will govern the spectrum. These charging effects will 
manifest themselves in jumps in the compressibility as a function of charge density. The 
simple dot model of [28] was able to account for this effect and even demonstrated that 
charging events will take place along lines of integer filling factor, in agreement with our 
calculations. An estimation for the cross-over from linear to non-linear screening can 
be found by very general considerations [70]. An insulating island where n(r) = uq is 
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Figure 10. Cross-sections of the system of Figure 9 at x = lOOnm and filling factors 
v = 0.1, u = 0.5, and v = 0.9 (bottom to top). Full curves correspond to n(r), broken 
lines show (22). A thick horizontal dotted line shows the average charge density, vno. 
The grey areas, bounded by thin dotted lines, indicate the complete band. 

confined by the force of the impurity potential, VVi(r), around its edge. Thereby, the 
Coulomb interactions opposing this force making the edge of the full region metallic. 
The size of the edges, i.e. the size of the metallic region, is then determined by the 
Coulomb force noe^/27reeo. Only if the Coulomb force acquires a magnitude comparable 
to the typical confining potential force, the metallic edges of the full islands will connect 
and dominate over the insulating regions. The typical force of our impurity potential is 
given by (|VVi(r)p) = ni{ws^)s/'^d'^- We would like to remark that with Ni = 288, the 
expected standard deviation of {ws^)s is ~ 2%, which makes the typical force a reliable 
characteristic of Vi for finite sample calculations. Prom equating the typical force with 
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Figure 11. Spatial distribution of HF-interacting electron density n{r)/no a.t B — 4T, 
7 — 0.3 and v ~ 0.9 as indicated by a color scale. Solid contour lines show the 
equipotential lines of (21). The broken lines indicate unscreenable (insulating) regions. 

the Coulomb force we can derive an expression for the minimal required density hq which 
corresponds to a minimum magnetic field -Bmin = noh/e below which linear screening 
breaks down for any density. Therefore, uq determines the width of the charging stripes 
An and we find 

Note that An is indeed independent of B and ne as observed in the experiments. 
In Figures 5, 6, and 7, we have indicated the breakdown of the linear screening 
regime by dashed white lines. The points at which the lines merge indicate Smin- 
Evidently, (23) nicely estimates the widths of the observed stripes for all three disorder 
configurations used. Furthermore, we have tested the criterion for breakdown by plotting 
compressibilities as a function of ne and disorder strength W/(P. Figure 12 shows the 
result for W/d"^ between ImeV and 3meV at B = 3.5T. The dashed white lines again 
indicate (23). In order to confirm the dependence of the stripes on the ratio between 
W and (i^ only, the plot has been divided into two regimes. Between ImeV and 2meV, 
we kept d = 40nm as a constant and varied W, and between 2meV and 3meV we kept 
W/nm"^ = 3.2meV constant whilst varying d, accordingly. The results confirm (23). 
Deviations from the expected behaviour can be explained with the proximity to the 



Compressibility stripes for mesoscopic QH samples 



19 



disorder dominated regime for higher values of W/d"^ where B ^ -Bmin- This regime is 
strongly disorder dominated and charging effects become much less pronounced at the 
band edges. 




1.0 1.5 2.0 2.5 3.0 



W/drimeVl 

Figure 12. Inverse electronic compressibility k^^ for a HF-interacting system 
(7 = 0.3) of size L = 300nm in the {W/(P,nc)-p\a,ne. Dashed white lines indicate 
expected boundaries of linear screening as calculated from (23), dashed black lines 
show integer filling factors. 

10. Conclusions 

We have investigated numerically how electron-electron interactions affect the 
localization properties of a 2DES under influence of a strong perpendicular magnetic 
field. We therefore diagonalized the Hamiltonian in a suitable basis and treated 
interactions as an effective mean field. 

Our calculations reveal substantial differences in the electronic compressibility 
between non-interacting and interacting systems when viewed as a function of magnetic 
field and carrier density. For interacting systems, we find strongly incompressible stripes 
of constant width around integer filling factors. We show the dependence of the width of 
the stripes on the disorder configuration and compute the width based on a force balance 
argument. These results are in very good agreement with recent imaging experiments. 
Moreover, we find strong (/-factor enhancement as well as negative compressibility in the 
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band centres, also consistent with experiments. We demonstrate that the incompressible 
patterns can be attributed to non-linear screening effects in the tails of the Landau 
bands. For magnetic fields larger than i?min, the effects of linear screening — and hence 
interactions — dominate in the {B, nc)-plane. Thus, our results support the existence of 
a greater variety of transport regimes due to electron-electron interactions in the integer 
quantum Hall effect. 

Similar compressibility patterns have also been observed around fractional filling 
factors u — 1/3, 2/5, and 2/3 [17]. Energy gaps at fractional filling, e.g. u — p/{2p+ 1), 
with p being an integer, are a consequence of electron correlations which are absent in HF 
approximation. However, with the formal analogy [71] between IQHE and FQHE put 
forward by the composite fermion (CF) model [72,73], let us venture a few statements 
about compressibility patterns around those fractional fillings. It is argued that the 
FQHE can be regarded as a manifestation of the IQHE for CFs in an effective magnetic 
field B* = Bi, — Bj,=i/2 [74, 75]. If we pretend to have obtained results for CFs in an 
{B*,ncF) plane, a transformation back to electrons yields an increase in the density of 
charging lines (per rie) by a factor of 2p-|- 1. Indeed, in the above mentioned experiment 
an increase of 3 has been found for i/ = 1/3. Furthermore, such a transformation predicts 
a dependence of the width of the incompressible stripes on the filling factor as well as a 
strong increase of -Bmin when fractional filling factors approach u — 1/2. This remains 
yet to be explored. 
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Appendix A. Plane wave matrix elements and boundary conditions 

For n >m, the plane wave matrix elements Sn,k;n',k'{Q) read 



(n^l e.m- rl\mj) = 5;,.-,y|^e-*+i-(^+^) (^-^^J'"^ L^^^ (|) , (A.l) 

where L^{x) is the generalized Laguerre polynomial. For periodic boundaries, the delta 
function is defined with modulus, i.e. 

, _ . 1 if mod(a - b, N^) = 0, 
^'^'-^ otherwise. ^^"^^ 
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The periodic boundary conditions require careful treatment of the overlap between 
Landau functions at opposite sites of the sample. We only take \mj) to be a replicated 
Landau function in one unit cell to the right and one to the left of the sample, whereas 
{nk\ remains in the base cell. A check of the implementation can be carried out by 
noting that a shift of the impurities by L/A^^ along the x or the y direction should only 
shift the wave-functions, ^a(r), by the same value. Finally we would like to remark 
that the complexity of summations involving plane wave matrix elements can be greatly 
reduced by neglecting terms where exp(— g^/2) < e. One can apply the restriction 
\qx\'^ < max(0, — 2 ln(e) — Qy) to summations over q^., usually leading to a substantial 
reduction of complexity even for e ~ lO"^'^. 

Appendix B. Bielectronic integrals 

Summations containing the bielectronic integrals 

Gn'w ^^v{q)Sn,k;n',k'{q)Sm,l;m',l'{-q) (B.l) 

can be substantially simphfied by virtue of the delta function (A. 2) contained in the 
plane wave matrix elements. Ultimately, two of the summations in (8) drop out. We 
can replace qyhjk — k' and I' by I + k — k', and if n > m as well as n' >m', we get 

, (^y'-"'(zl^)--^.-.(£),_(p. ,3.3) 
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